Spin Helix of Magnetic Impurities in Two-dimensional Helical Metal 
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We analyze the Ruderman-Kettel-Kasuya-Yosida(RKKY) interaction between magnetic impuri- 
ties embedded in the helical metal on the surface of three-dimensional topological insulators. Apart 
from the conventional RKKY terms, the spin-momentum locking of conduction electrons also leads 
to a significant Dzyaloshinskii-Moriya (DM) interaction between impurity spins. For a chain of 
magnetic impurities, the DM term can result in single-handed spin helix on the surface. The hand- 
edness of spin helix is locked with the sign of Fermi velocity of the emergent Dirac fermions on the 
surface. We also show the polarization of impurity spins can be controlled via electric voltage for 
dilute magnetic impurity concentration. 



PACS numbers: 72.25.Dc, 73.20.-r, 75.30.Hx, 85.75.-d 

It is now known that a three-dimensional (3D) insula- 
tors with time reversal symmetry can be classified into 
two categories: ordinary and topological insulator(TI) 
[IH5] . Though both of them are fully gapped in the bulk, 
TI has metallic surface states robust against any time 
reversal invariant perturbations. These metallic surface 
states dubbed as "helical metal" can be described by a 
two-dimensional(2D) massless Dirac equation. Compar- 
ing to other emergent relativistic materials such as single- 
layer graphene, the surface states of TIs have two unique 
features. One is that the number of Dirac nodes is odd. 
In fact, according to a no-go theorem|S], the Dirac nodes 
always appear in pairs in the conventional 2D lattice such 
as graphene. A surface state of TI can invalidate no- 
go theorem because a pair of Dirac nodes are separated 
onto two opposite surfaces. TI with single Dirac cone 
has been theoretically studied and experimentally ob- 
served by angle-resolved photo-emission spectroscopy |7l- 
HTj . The other feature is that electron spin and mo- 
mentum are intimately locked in helical metal, origi- 
nated from strong spin-orbit coupling. Evidence of the 
spin-momentum locking has also been observed in recent 
experiments [12H15] . 

The effect of magnetic impurities in a metal is a cen- 
tral issue in condensed matter physics, since the cou- 
pling between impurity and conduction electrons can 
lead to many intriguing phenomena such as Kondo ef- 
fect and Ruderman-Kettel-Kasuya-Yosida(RKKY) inter- 
action. The analogy of these effect in a helical metal of- 
ten reveals intrinsic property of a TI, which has received 
considerable attention recently [ToT - fH)] . In this letter we 
point out a novel feature in RKKY interaction between 
impurity spins in helical metal, which manifests directly 
the two features of TI mentioned above. Firstly, due 
to the spin-momentum locking of conduction electron, 
we show the induced RKKY interaction must contain 
an anisotropic Dzyaloshinskii-Moriya(DM) term by both 



symmetry analysis and perturbation calculation, which 
is usually absent in systems like graphene [201 - F22] . Since 
the emergent Dirac Hamiltonian of helical metal is a pure 
spin-orbit coupling, the magnitude of DM term is of the 
same order of the conventional RKKY one. Secondly, 
the DM term can lead to single-handed spin helix if the 
impurities are aligned into a chain, of which the handed- 
ness depends on the sign of the Fermi velocity vf of the 
helical metal. For strong TIs with a single Dirac cone 
attached to each surface, vp is fixed, therefore the cor- 
responding spin helix must be single-handed. The spin 
helix is illustrated in Fig. [I] where two types of Dirac 
fermions are considered and the spin helix are perpen- 
dicular(Fig. [iji) and parallel (Fig.[T]3) to impurity chain, 
respectively. This is a hallmark of a strong TI, and could 
be possibly detected by spin-polarized scan tunneling mi- 
croscopy, or optical measurement like Kerr effect. 




FIG. 1. Illustration of spin helix of a chain of magnetic im- 
purities. Two types of Dirac Hamiltonian, H§ and //o(see 
Model section), are considered, which lead to spin helix (ma- 
genta vectors) perpendicular(a) and parallel(b) to the chain 
direction, respectively. The insets show spin orientation of 
conduction electron in momentum space. 

Model: We start from the Hamiltonian of the conduc- 
tion electrons with only one Dirac cone 

6 s = 4(^ k • " - E F i)c k , (1) 
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with Fermi energy Ep and Pauli matrix er. Via a unitary 
transformation aks = [ckt + sgn(i;p')se _i6 ' k Ck^]/v / 2 with 
6*k the angle of vector k and s = ±1, it can be diagonal- 
ized as Hq = J2 s =±i( s ^\ v f\^ ~ ^/)^Li^ks- In practice, 
the recently discovered TI material, e.g., Bi 2 Se3 has sur- 
face state described by the Rashba type of Hamiltonian 
= J2k c k [^F(k x cr) z — £^/]ck[S E], which is related 
to Hq by a rotation of k around z axis by it/2, hence we 
can focus on Hq in the following. 

Suppose there are N imp impurity spins located at r n 
denoted by S„, which interact with conduction electrons 
via the following spin-spin coupling 

N imp 

Hi=J2 A z ^(r„)^ + A±[r(r n )^ + ^(r„)^],(2) 

n=l 

where s x ' y ' z (r) are the spin operators of the conduction 
electron at r. The coupling constants are assumed to be 
isotropic in xy plane A^ = X y = \±, but may differ from 
the z-component. The total Hamiltonian is simply the 
sum H = H„ b + Hj. 

Before proceeding with the discussion of the effective 
interaction between impurity spins, we follow RKKY[2"3T- 
[25] to divide Hj into two parts: Hj and Hji. Hj con- 
sists of the diagonal term of a ks ak s and can be written 

as H I0 = A± J]k( e k ' «k) (ek ■ S imp ), where e k = k/k, 

s k = c k o-c k , and S imp = V -1 Yln=7 S n is the density 
of impurity spin. And the remaining part is off-diagonal 
denoted by H\\. Notice that for the helical Hamiltonian, 
the electric current of the conduction electrons is propor- 
tional to their spins, for instance, J = evF X)k ®k for Hq , 
therefore Hj actually implies a direct coupling between 
the electric current and the impurity spins, which pro- 
vides a mechanism to control the impurity spins through 
the electric voltage for dilute impurity concentration. 

RKKY interaction: For simplicity we focus on the case 
of two impurities at ri and r 2 , respectively. The results 
given below via the symmetry analysis are also valid for 
many impurity case as long as we consider two-body in- 
teractions. The RKKY interaction can be obtained by 
integrating out the degree of freedom of conduction elec- 
trons, followed by an expansion to the second order of the 
coupling parameters A Zi ±[;23 25 . The effective interac- 
tion between impurity spins thus obtained contains only 



the even order terms for time reversal invariant system. 

To the lowest order, the most general form of the in- 
teraction is H rkky (v 12 ) = Ylaiaa F "i"2 i^S" 1 S 2 2 with 
r 12 = ri — r 2 . Without spin-orbit coupling as in the usual 
case, only terms with a\ — a 2 can exist due to a rota- 
tional symmetry in spin space. However, the situation is 
quite different for the helical metal, since Hq is invari- 
ant only under a joint rotation in both spin and orbital 
space, which allows more terms as soon become clear. 
Let us consider a global rotation around z axis by <p with 
the form R e ,, v = exp[ipe z ■ s cond ] x exp{itpe z ■ (Si + 
S 2 )] x exp[iipe z ■ L] with total spin of conduction elec- 
trons s CO nd = J drc 1 (r)crc(r) /2 and orbital angular mo- 
mentum L = J dr&(r)(r x p)c(r). Under this rotation, 
Hq is obviously invariant, but Hj(ri,r 2 ) — > Hi(r' l ,r' 2 ) 
with r' = e -*v(»Pv-vp») re i<p(a*y-»p*) ) w hich breaks the 
rotational symmetry. However since the energy behaves 
like a scalar, one expects H rkky (ri 2 ) = H rkky (r' 12 ), then 
what we are left with is to construct R e , v -invariants with 
vectors e 12 (= i"i 2 /ri 2 ), Si and S 2 . 

There are many such R ez ^-invariants terms, most of 
which can be further eliminated by the following two 
symmetries. The first one is exchanging ri and r 2 which 
obviously leaves H invariant. Thus terms like e 2 ■ (Si x S 2 ) 
are not allowed. The second one is a global rotation 
R ei2l 7r around axis ei 2 by n, which rules out invariants 
like [ei2 x (Si x S 2 )] z that changes sign under rotation 
Re 12 ,7r- Hence, we are left with the following two terms 
in addition to the conventional ones, ei 2 • (Si x S 2 ) and 
(ei2 • Si)(ei2 • S 2 ). One then constructs the effective 
Hamiltonian of impurity spins as 

H rkky (r 12 ) = r (ri 2 )[A|(^ + SfS v 2 ) + X 2 Z S^} 
+ri(ri 2 )A|[(ei 2 -Si)(ei 2 -S 2 )] 
+r D M(ri2)A 2 A±[ei 2 • (Si x S 2 )], (3) 

where the To-term is just the conventional RKKY terms. 

Derivation of Tq.i.dm'- The coefficients Tq^jjm in 
Eq. Q can be computed following RKKY's second or- 
der perturbation [23 25 . Let us consider the following 
generalized partition function which is partially traced 
over conduction electrons, Z — Tr conc ie~ ,3 (- Ho+ - Hl ). To 
the quadratic order of A Zl ±, we obtain 



|^kEEE V (1 / k2S2) t E ^Ks^sM E A b < S2;klsl ^]e^--)^-^)} (4) 

[ y k lSl k 2 s 2 n,m ?k i s i a=x,y.z b=x,y,z J 



where Zq is the partition function of conduction elec- tion function, and the three -F-functions take the forms 
trons, / ks = [1 + e /3(*ft|«F|k-%)]-i i s the Fermi distribu- ^k lSi; k 2S2 = \sgn(v F )( Sl e~ iek i + s 2 e l9 ^) i F^.^ = 
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lsgn(v F )( Sl e-^ - s 2 e ie *2) and F£ iaiika , 3 = |(1 - 
SiS2e~^ 6 ' k i _9k 2)), which depend only on 9^ 2 , the angle 
of ki^. By comparing with Eq. ([3]), one can extract the 
coefficients Fo,i,d m fromEq. Q straightforwardly. After 
integrating over 9k 12 , we have T (r) = —2[A(r) — B(r)}, 
Ti(r) = -4B{r) and T DM (r) = -2C(r), where the func- 
tions A, B and C read in the zero temperature limit 

A ^ = 8h\v F \r 3 ^ J x i dx i J X2dx 2 3o(xi)So(x 2 ) 
9(xf - xi)9(x 2 - xf) 

S 2 X 2 - SiXi 

B(r) = 1 [ xxdxx I x 2 dx 2 di(xi)di(x 2 ) 

8H\v F \r 3 f^J J 

9(xf — xi)9(x 2 — xf) 

S 2 X 2 - SiXi 
C ^ = S^F$^j XldXl j X2dX2 

[sidi{xi)do(x 2 ) - s 2 3o(xi)di(x 2 )] 

9{x F - x 1 )9(x 2 - x F ) 

x (5) 

s 2 x 2 - si^i 

where x = kr, xp = kpr, 9(x) is Heaviside step func- 
tion, and 3o,i( x ) are the first and second order Bessel 
functions, respectively. Note that r (r) has already been 
obtained in Ref.|18j. which agrees with ours. It is the Ti 
and T dm terms that have not been addressed before, and 
the physics discussed in this letter is essentially coming 
from them. 

The numerical results of T Q iDM( r ) are plotted in 
Fig. [3]where a cutoff is set for x and the Cauchy principle 
integration is understood when singularities are encoun- 
tered in Eq. ([5| . Despite of the oscillation of the conven- 
tional RKKY interaction with respect to r for large kp, 
it is interesting to note that the DM interaction is of the 
same order of the conventional ones, and depends on the 
sign of vf as seen from the expression of C(r) in Eq. 

Single-handed Spin Helix: To reveal the effect of DM 
interaction, we focus on a simple case of one-dimensional 
chain of impurities aligned with equal spacing d along 
x-axis, and the corresponding Hamiltonian reads 

Hi m p = y ^ (—J z S n S n+1 — JyS^Sf l+1 + J x S n S n _^_i) 
n 

+ JDM(SnSf l+1 — S^S Z n+l ) (6) 

where J z = -T (d)X 2 z , J y = -T (d)X 2 ± , J x = [T {d) + 
Ti(d)]\±, and Jdm = —^DM(d)X z X±. Here we consider 
the case for kp close to the Dirac point, where we can 
simply assume J y>z > J x > as can be read off in Fig. [2Jd 
for r 1. The other situations including two dimensional 
arrays of magnetic impurities are left for future studies. 

In this case, the spins are lying in yz plane in the 
classical limit. Therefore, we can introduce two variables 
9 n and p n to describe the spin rotation in j/z-plane and 
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FIG. 2. Plot of To,i,dm as functions of the distance between 
impurities r, /cfOo = 1-8 (left) and /cfcio = 0.5 (right). The 
distance between impurities is in unit of lattice spacing oq. 
The coefficient Shvpag is taken as one for convenience. 

the amplitude fluctuation along x-axis, which satisfy the 
commutation relation [9 ni p m ] — iS m . n , and are related to 
S„ through S* = p, S% + iS* = e l§ [{S~p){S + p+l)Y' 2 
and SI - iS z n = [(S - p)(S + p+ l)] 1 ^-^. One may 
verify these equations satisfying the angular momentum 
algebra, and S 2 = S(S +1). In continuum limit the 
effective Hamiltonian for the dynamics of impurity spins 
can be written as 

H imp = J dx[-D 1 cos(ddj) + D 2 cos(20) 

+D 3 sm(ddj) + D p 2 } (7) 

with [9{x), p(x')\ = iS(x — x'). The parameters are given 
by D = J x d, Di = (J z + J y )S(S + l)/(2d), D 2 = (J z - 
J y )S(S + l)/(2d) and D 3 = J DM S{S + l)/d. 

We first consider the isotropic case, i.e., J z — J y - The 
9(x) field has a helical background(see Fig. [lja)), with 
helical angle r\ satisfying tan(?yd) = 2 Jdm /(Jz + Jy)- The 
sign of r\ is the same as that of Jdm, and its amplitude 
is saturated to it /2d as Jdm is large enough as shown in 
the lower panel in Fig. k5[ Replacing 9{x) in Eq. with 
rjx+9(x), we have H imp = J dx[—^D\ + D 2 cos(dd x 9) + 
Dq/} 2 } describing the quantum fluctuation upon the he- 
lical background. The low energy excitations can be 
obtained by expanding cos(dd x 9) to the second order 
of d, which leads to a linear spectrum in k as ujk — 
d{2D fl 2 {D 2 + D 2 Y/ i \k\. 

Next we consider the effect of spin anisotropy, i.e. 
Jz 7^ Jy Without loss of generality, we assume J z > J y (if 
J z < J y , one can shift 9 — > 9 + n/2 to get posi- 
tive D 2 ). In this case the Hamiltonian becomes a sine- 
Gordon(SG) model H mp « / dx^[K- 1 (d x 9) 2 + Kp 2 ] + 
hd x 9 + D 2 cos(y/8Tr9), which is written in the standard 
form by rescaling 9 and p in order to use the well known 
results in literature, e.g. in Refs.|261 \Tl\ . The coeffi- 
cients take the form u = d^J J x (J y + J Z )S(S + 1), K = 
■K~ 1 y/J x /(J Z + J y )S(S + l) and h = V^JdmS(S + 1). 
In case K > 1, the cosine term is irrelevant, so that the 
system is massless and the spin helix takes place for any 
finite h. However in case K < 1, the theory turns out 
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FIG. 3. A schematic plot of phase diagram of J z — J y vs. 
J dm in the upper panel for K < 1, and helical angle as a 
function of Jdm in lower panel. The left helix, non-helical 
state and right helix regions belong to different topological 
sectors of sine-Gordon model with negative, zero, and positive 
topological charges, respectively. 



to be massive, and a critical value of |Jdm| exists, only 
above which the spin helix can occur in forms of mas- 
sive soliton excitations of 8 field [26] [27] which connect 
different classical vacua of the SG model. In our case, 
since < J x < J y , z , we are in the region of if < 1. 
A schematic phase diagram for this case is given in the 
upper panel of Fig. [3j where the left helix, non-helical 
and right helix regions correspond to different topologi- 
cal sectors of sine-Gordon model with negative, zero, and 
positive topological charges, respectively. 

Our analysis of DM interaction so far is focused on 
Hq, where the spin helix is in yz plane perpendicular 
to the impurity chain. For helical metal described by 
H$, the RKKY interaction can be obtained by rotat- 
ing ei2 around z-axis by 7r/2 in Eq. ([3]). As a con- 
sequence, the r -term is invariant, r^term becomes 
ri(ri 2 )A|(ei2 x Si) 2 (ei2 x S 2 ) z and T^M-term becomes 
^DM{ri2)^±^z[ei2 x (Si x S 2 )] 2 . Notice that (e 12 x 
Si) 2 (e 12 x S 2 ) z = SfS% + SfS y 2 - (e 12 • S^e^ • S 2 ), 
therefore Ti-term simply changes the coefficients of the 
first two terms of Eq. (|3| . Only the change of Yd M-term 
is essential, which makes the spin rotate in xz plane in- 
stead of yz plane as illustrated in Fig. [T|d. 

Summary: The surface state of TI is metallic with 
strong spin-orbit coupling, in which the magnetic impu- 
rities coupled with conductance electrons can be polar- 
ized by the electric voltage. There is also an effective 
interaction between impurity spins mediated by the con- 
duction electrons, which includes a DM interaction with 
the same order of amplitude of the isotropic RKKY in- 
teractions. For ID chain of impurities, this could lead to 
a single-handed spin helix, and the handedness is locked 
with the sign of Fermi velocity vp of the emergent Dirac 



particles. 
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